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The mixing between sterile and active neutrinos is taken into account in the calculation of Big 
Bang Nucleosynthesis. The abundances of primordial elements, like D, 3 He, 4 He and 7 Li, are 
calculated by including sterile neutrinos, and by using finite chemical potentials. It is found that 
the resulting theoretical abundances are consistent with WMAP data on baryonic densities, and with 
limits of LSND on mixing angles, only if 7 Li is excluded from the statistical analysis of theoretical 
and experimental results. 
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I. INTRODUCTION 



In recently published papers [j], 0| the sensitivity of the 4 He primordial abundance, upon distortions of the light 
i - 1, neutrino spectrum induced by couplings with a sterile neutrino, was analyzed. The effects due to the mixing between 
*-£h ■ sterile and active neutrinos [H, HJ reflect upon Big Bang Nucleosynthesis (BBN) in a noticeable manner. Previous 
studies on this matter can be found in 0, 0, @, where the effects of mixing upon BBN in presence of primordial 
leptonic asymmetry have been investigated, and stringent limits on the mixing due to BBN have been presented. 
' Similar studies have been presented in [f| @]. A review on inclusion of sterile neutrino in cosmology was presented in 
!• _ 

i The results of [2| may be taken as a solid starting point for a systematic analysis of the sterile-active neutrino 

mixing upon cosmological observables, like the primordial abundance of light elements. By the other hand, the 
mixing mechanism between sterile and active neutrinos has been studied in detail (see Q and references therein), so 
that the calculation of neutrino distribution functions can readily be performed. The information about the neutrino 
distribution function, in the flavor basis and at a given temperature, is an essential element in the calculation of the 
neutron decay rate, which is a critical quantity entering BBN (if! [Tl|. 
. Direct physical consequences upon BBN, due to the mixing between active and sterile neutrinos, have been explored 
in [13, [H| . Following the arguments presented in [12] , and in the framework of the standard cosmological model, 
sterile neutrinos would produce a faster expansion rate for the Universe and a higher yield of 4 He. This is, indeed, a 
severe constraint on neutrino mixing since a higher predicted abundance of 4 He may be in conflict with observational 
data Another constraint on active-sterile neutrino mixing is the neutrino mass derived from Cosmic Microwave 
Background Anisotropy (CMB) [T3|.The analysis of constraints presented in [l5| focus on the mixing scheme at the 
level of the neutrino mass hierarchy, and it suggests the adequacy of the non-degenerate mass hierarchy to set limits 
on the mass difference between active and sterile neutrinos, 5m\_ s . The study of [l5| confirms the notion about the 
convenience of the three active + one sterile neutrinos scheme. 

In standard BBN calculations, the mixing of sterile and active neutrinos affects the leptonic fractional occupancies, 
which are essential quantities appearing in the expression of the weak decay rates. Thus, one needs to know, as input 
of the calculations, the parameters of the proposed mixing scheme, the neutrino mass hierarchy and the leptonic 
densities [l3| . With these elements one can calculate neutron-decay-rates and neutron abundances, by assuming the 
freeze-out of weak interactions [lfj, |12J. The effective number of neutrino generations, N„, is fixed by the analysis of 
CMB 0, [l?], HI]. Current limits on the neutrino degeneracy parameter, for light (electron) neutrinos, r\i runs 
from —0.1 to 0.3 [l7l. Il8j|. For a detailed presentation of the formalism, in the context of relic-neutrino asymmetry 
evolution see 

In this work we focus on the calculation of the abundances of D, 3 He, 4 He and 7 Li, in presence of sterile-active 
neutrino mixing in the three flavor scenario, and for the normal and inverse neutrino mass hierarchies Q. We have 
compared the calculated values with data [lj| [2(| [H| and determined the compatibility between them by performing a 
\ 2 statistical analysis. Since the theoretical expressions depend on the mixing angle sin 2 2(f>, the square mass difference 
<5mf 4 (normal mass hierarchy) or 5m \ A ( inverse mass hierarchy), and the baryonic density fish 2 , we have adopted 
the LSND limits on the mixing angle [22l. l23l. l24j and the WMAP results on the baryonic density [2||, as constraints. 

The paper is organized as follows. In Section [II] we briefly present the essentials of the formalism. Section HTT1 is 
devoted to the calculation of the neutron decay rate and BBN abundances. In Section [IV] we present and discuss the 
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results of the calculations. Conclusions are drawn in Section fVl 



II. FORMALISM 



The mixing between active neutrino mass eigenstates i>i (i = 1,2,3), leading to neutrinos of a given flavor i>k 
(k = light, medium, heavy), is described by the mixing matrix U (26j 
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where cy and s 7 ; 



stand for cos 0y and sin0y, respectively, and CP conservation is assumed [261 ] - To this mixing we 
add the mixing of a sterile neutrino with: a) the neutrino mass eigenstate of lowest mass in the normal mass hierarchy, 
Vi, and b) to the one of the inverse mass hierarchy, z/3, by defining the mixing angle cj>, such that the new mixing 
matrix U is redefined as U ((f)) Q 
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for the normal mass hierarchy, and, 
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for the inverse mass hierarchy. The mixing between neutrino mass eigenstates, and particularly the inclusion of the 
sterile neutrino as a partner of the light neutrino, affects the statistical occupation factors of neutrinos of a given flavor. 
The equation which determines the structure of the neutrino occupation factors, in the basis of mass eigenstates and 
for an expanding Universe, can be written [27| : 
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where t is time, H is the expansion rate of the Universe, defined as H = y' 45 j^ — T 2 = fipT 2 , T is the temperature, 

E v is the energy of the neutrino, and Hq is the unperturbed mass term of the neutrino's Hamiltonian in the rest 
frame. The initial condition is fixed by defining the occupation numbers at the temperature Tq = 3 MeV [28| . 
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for the inverse mass hierarchy. 

To obtain the solutions of Eq.JH) we have written the source term, that is the commutator in the r.h.s of Eq.(|3]), in 
terms of the square mass differences, 6m 2 : = mf 
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The value of the mixing angle 6*13 is constrained by the upper limit given by [2!| so that tan #13 < 10 3 . The solution 
in the basis of mass eigenstates is 
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where the normalization constants are fixed by the initial conditions (T = Tq). The formal solution for the occupation 
number in the flavor basis, for the light neutrino flavor in the normal mass hierarchy, is: 
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in the inverse mass hierarchy. 

In the above expressions, rj is the ratio between the neutrino chemical potential and the temperature. This parameter 
depends on the adopted value of the leptonic number L [J Q • Explicit expressions of r\ versus L can be found in [l[ . 
In the present context we have taken 77 as an input for the calculations (see section [TIT)) . 



III. DECAY RATES AND NEUTRON ABUNDANCE 



In the following we shall outline the main steps of the calculation of neutron decay rates, for the electroweak 



processes n + e H 



p + v 



and 



p + e 



The starting point is the calculation of the reduced rates A± [10| 
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and the total neutron to proton decay rate 
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at lowest order in the quantity e v . For the sake of convenience we have introduced the more compact notation 
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(E 1 ^, T)| with a = c 2 3 c 2 2 for the normal mass hierarchy and a — s 2 3 for 



the inverse mass hierarchy, respectively. The function g(E v ,T) is the factor which contains the temperature T and 

The details of the calculations have 
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the energy E v in Eqs.fll]) and pH|) . and the variable y is defined as y = 

been discussed elsewhere [3(|, for the case of two neutrino mass eigenstates. The final expression for the neutron to 
proton decay rate is obtained by fixing the normalization Ao, of Eq. (|13[) . from the neutron half- life 
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in units of sec . Following Ref.[10(, the neutron abundance, until the freeze-out of weak interactions, is expressed in 
terms of the neutron to proton decay rate, X np of Eq. (|15[) as 

^neutrons - J^dW e™+" ( l + ^ ) XT ") W«) (lg) 

The quantity A nGU t r0 ns is, therefore, a function of X np and, consequently, of the occupation factors which contain 
the information about the mixing between active and sterile neutrinos. The next step consists on the calculation 
of primordial nuclear abundances. The method to calculate the BBN abundances was presented in Ref [ll|. It is a 
semi- analytic approach based on the balance between production and destruction of a given nuclear element, which 
requires the knowledge of A noutr ons- For details we refer the reader to [PH ]. 

The above presented framework shows that the calculation of primordial abundances may indeed be taken as a 
tool to test leptonic mechanisms, like the mixing between sterile and active neutrinos, as it has been pointed out by 
Kishimoto et al. 0. 

IV. RESULTS AND DISCUSSION 

To perform the calculations we have adopted the oscillation parameters determined from SNO, SK and CHOOZ 
measurements 29] . The mixing with the sterile neutrino, represented by the mixing angle </>, is taken as an unknown 
variable, within the limits fixed by the LSND data [22|, l23l. |24j. The mass splitting 8m\ A (or 5m| 4 ) was taken from 
the analysis given by Keranen et al. 0. The actual value is fixed at 5m? = 10~ n eV 2 . We have then calculated the 
neutron abundance, by applying the formalism of the previous section. The primordial abundances of D, 3 He, 4 He 
and 7 Li, have been calculated as described in [llj]. The baryonic density fish 2 (see Ref. [30]) was varied within the 
limits 0.010 < Vlsh 2 < 0.035. Concerning the value of r] we have varied it in the interval determined by the allowed 
values of the potential lepton number, C = 2L Ve + L v + L„ T , that is 0.0 < C < 0.4 [HQ. In the present calculations 
we have adopted the values 0.0 < r] < 0.07 which are consistent with the densities 0.0 < L„ c < 0.05 [l|. 

To determine the allowed values of the mixing angle </> we have performed a ^-minimization, after computing the 
primordial abundances. The data have been taken from Refs.[l9l. l2Ct l2l|. The results are shown in Figure 1, insets 
l.(a) and l.(b). The curves are the contour plots for results with comparable values of x 2 ■ Figure l.(a) shows the 
results obtained by the x 2 -analysis of theoretical and experimental values [HI H(J HH , including data on 7 Li. Figure 
l.(b) shows the results of the statistical analysis performed with the exclusion of 7 Li. In the first case, Figure l.(a), 
the absolute minimum is located at sin 2 2(f) = 0.000 ± 0.026, and Qsh 2 — 0.0253 ± 0.0015, both set of results have 
been obtained by using the solution ([9]) for the occupations. The smallness of the mixing angle does not contradict 
LSND results [22, HH, H3] , but the value of the baryonic density is outside the limits determined by WMAP [25| , that 
is: (^b^ 2 )wmap = 0.0223 ±0.0008. This disagreement between theory and data may be caused by large uncertainties 
in the 7 Li-data. As pointed out by Richard et al. [3l|, the validity of the data on 7 Li may be questioned by the 
uncertainties inherent to the physics of 7 Li in the interior of the stars, i.e; the turbulent transport in the radiative zone 
of stars. In contrast, the situation improves if the data on 7 Li are removed at the time of performing the statistical 
analysis (Figure l.(b)). For this case, the best value of the mixing angle is sin 2 2<p = 0.018 ± 0.098, and the baryonic 
density corresponding to the minimum, ilp/i 2 = 0.0216 ± 0.0017, is indeed consistent with the WMAP data. The 
anomalous feature associated with the inclusion of 7 Li in the set of data persists if other elements are removed from 
the data. We have verified it by systematically removing, one at the time, the abundances of D, 3 He, and 4 He, and 
keeping the data on 7 Li. In all cases the location of the minimum lies closer to the one of Figure l.(a). For the case 
of inverse mass hierarchy the occupation factor (|L0[) is strongly constrained by the value of #13 and the difference 
with respect to the thermal occupation factor vanishes. It means that the contour plot shows, for the inverse mass 
hierarchy, parallel lines to the vertical axis (sin 2 2<j>), since all possible values of sin 2 2<f> are allowed by Eq. (fTUj) when 
sin 2 — > 0. For both the normal and inverse hierarchy solutions, Eq. Q and (|10[) . particle number conservation 
was enforced, on the average, by the factor /ip (see its definition following Eq.(|3])). Because of the high temperature 
we have not included collision terms in Eq. (]3|). 

Similar results, related to the abundance of 7 Li, have been obtained in the calculations of nuclear abundances in 
the context of cosmological models [32], [3^, [34[ , and also in the case of a two neutrino mixing [3(| ■ 

To investigate the dependence of the above results on the parameter 77, we show in Figure [2] the values of the mixing 
angle sin 2 2<f>, obtained from the x 2 analysis, as functions of the chemical potential. The calculations have been 
performed by excluding the data on 7 Li. Our present results are very much in agreement with the results reported 
in Ref. since the values shown in Figure [2] display a small variation in a relatively large domain of values of r\. 
Finally, the best values of the baryonic density and the mixing angle, both with and without including 7 Li in the 
analysis, are shown in Table [J as functions of r\. In agreement with the expectations of [To| . and with our owns, both 
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FIG. 1: Statistical analysis of the calculated nuclear abundances. The curves are the contour plots for results with comparable 
X 2 values. The contours correspond to increasing values of x 2 , from bottom to top. The calculations were performed by taking 
the sterile-active neutrino mixing, sin 2 2<j>, and the baryonic density, Qsh 2 , as variables. The inset l.(a) shows the results 
obtained by the \ 2 analysis of theoretical and experimental values, including data on 7 Li. Inset l.(b) shows the results of the 
statistical analysis performed with the exclusion of 7 Li. The results shown in the figure have been obtained with the solution 
corresponding to the normal mass hierarchy. 



sets of results do not differ much, or at least they do not show a pronounced dependence, with respect to the chemical 
potential. 





All data 
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sin 2 2<f) 


0.00 


0.0253 ± 0.0015 


0.000 ± 0.026 


0.0216 ±0.0017 


0.018 ±0.098 


0.01 


0.0250 ± 0.0014 


0.000 ±0.010 


0.0216 ±0.0020 


0.002 ±0.022 


0.02 


0.0248 ± 0.0014 


0.000 ± 0.015 


0.0218 ±0.0020 


0.004 ± 0.030 


0.03 


0.0246 ± 0.0012 


0.000 ± 0.034 


0.0216 ±0.0018 


0.008 ± 0.080 


0.04 


0.0244 ± 0.0016 


0.000 ± 0.039 


0.0216 ±0.0019 


0.018 ±0.090 


0.05 


0.0244 ± 0.0016 


0.000 ± 0.056 


0.0216 ±0.0017 


0.052 ±0.090 


0.06 


0.0244 ± 0.0016 


0.000 ± 0.101 


0.0216 ±0.0018 


0.108 ±0.090 



TABLE I: Best values of the mixing angle and of the baryonic density, determined from the \ analysis of the calculated 
abundances, as functions of the parameter rj. Left and right sides of the table show the results obtained with and without 
considering the data on 7 Li, respectively. 
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FIG. 2: Best values of the mixing angle, sin 2 2(j>, determined from the \ 2 analysis of the calculated abundances, as functions of 
the parameter r\ 



V. CONCLUSIONS 

In this work we have calculated BBN abundances, by including the mixing between active and sterile neutrinos. As 
pointed out by Kishimoto et al, the BBN abundances are sensitive to active-sterile neutrino mixing, indeed. Kishimoto 
et al. Q have demonstrated the sensitivity of the 4 He abundance on the distortion of the light neutrino spectrum 
produced by the mixing with a sterile neutrino. In our case, the statistical analysis of the compatibility between 
theoretical and observed nuclear abundances, indicates the existence of a clear sensitivity of the results upon active- 
sterile neutrino mixing, too. In performing our analysis, we have considered the WMAP baryonic density, together 
with the LSND constraint on the sterile-active neutrino mixing. The comparison between calculated and observed 
abundances indicates some sort of anomaly in the abundance of 7 Li. Similar difficulties, related to the determination 
of the abundance of 7 Li, have been reported previously (3lj |. in the context of the physics of the interior of stars. 
We found that the consideration of the abundance of 7 Li, in presence of active-sterile neutrino mixing, excludes the 
WMAP value of the baryonic density. This exclusion is not observed when the other nuclear abundances are not 
included in the analysis. 
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